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Abstract 

Background: Chemicals generically referred to as flavonoids belong to the group of phenolic compounds and 
constitute an important group of secondary metabolites due to their applications as well as their biochemical 
properties. Flavonoids, which share a common benzo-y-pyrone structure, constitute a kind of compound which are 
highly ubiquitous in the plant kingdom. 

Findings: The M06 family of density functionals has been assessed for the calculation of the molecular structure and 
properties of the Naringenin flavonoid. The chemical reactivity descriptors have been calculated through Conceptual 
DFT. The active sites for nucleophilic and electrophilic attacks have been chosen by relating them to the Fukui 
function indices and the dual descriptor f (2) (r). A comparison between the descriptors calculated through vertical 
energy values and those arising from the Koopmans' theorem approximation have been performed in order to check 
for the validity of the last procedure. 

Conclusions: The M06 family of density functionals (M06, M06L, M06-2X and M06-HF) used in the present work leads 
to the same qualitatively and quantitatively similar description of the chemistry and reactivity of the Naringenin 
molecule, yielding reasonable results. However, for the case of the M06-2X and M06-HF density functionals, which 
include a large portion of HF exchange, the calculations considering the validity of the Koopmans' theorem lead to 
negative electron affinities. 
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Findings 

Introduction 

Chemicals generically referred to as flavonoids belong 
to the group of phenolic compounds and constitute 
an important group of secondary metabolites due to 
their applications as well as their biochemical proper- 
ties. Flavonoids, which share a common benzo-y-pyrone 
structure, constitute a kind of compound which are highly 
ubiquitous in the plant kingdom. Over 4, 000 different 
naturally occurring flavonoids have been discovered, and 
only in the case of flavones, a specific type of flavonoids, 
over 36, 000 different chemical structures are possible. 
Flavonoids are present in a wide variety of edible plant 
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sources, such as fruits, vegetables, nuts, seeds, grains, tea 
and wine [1]. 

The knowledge of reactivity on a molecule is an essen- 
tial concept; it is of a crucial interest because it allows 
to understand interactions that are operating during a 
reaction mechanism. In particular electrostatic interac- 
tions have been successfully explained by the use of the 
molecular electrostatic potential [2,3]. 

On the other hand, there is no a unique tool to quan- 
tify and rationalize covalent interactions, however since 
2005 a descriptor of local reactivity whose name is simply 
dual descriptor [4,5], has allowed to rationalize reaction 
mechanisms in terms of overlapping nucleophilic regions 
with electrophilic regions in order to get a maximum sta- 
bilization thus leading to final products or intermediates; 
all those favorable nucleophilic-electrophilic interactions 
have been explained as a manifestation of the Principle 
of Maximum Hardness [6] in addition, chemical reactions 
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have been understood in terms of the The Hard and Soft 
Acids and Bases Principle [7-10], principle that has been 
used even with the aim of replacing the use of the Molec- 
ular Orbital Theory to understand the whole Chemistry 
[11]. In fact the present work is a good chance to test the 
capability of the most recent reactivity descriptors coming 
from the Conceptual DFT [12-15], therefore the frame- 
work of this conceptual theory will be presented in the 
next section. 

Naringenin (C15H1205; mol. wt. 272.3; IUPAC (2S)- 
5,7-Dihydroxy-2-(p-hydroxyphenyl)-4-chromanone, is a 
plant flavonoid that is extracted from citrus fruit. 
Its properties are estimated as follows: Polarizability = 
27.3945 i 3 , Molar Refractivity = 70.619 cm 3 /mol, Polar 
Surface Area = 86.99 i 2 , vdW volume = 234.6628 i 3 , 
logP = 3.398, and Complexity = 476.8368. It must be kept 
in mind that the objective of this work is not to perform 
an evaluation of the antioxidant properties of Naringenin, 
but to do a comparative study of the performance of 
the M06 family of density functionals for the description 
of the chemical reactivity of this prototypical flavonoid 
whose molecular structure is shown in Figures 1 and 2. 

Theory and computational details 

Morell et al. [5,9,11,16-19] have proposed a local reactiv- 
ity descriptor (LRD) which is called the dual descriptor 
(DD) f^ 2 \v) = A/(r). In spite of having been discovered 
several years ago, a solid physical interpretation was not 
provided in such a moment. [20]. They used the notation 
A/(r), but currently it has been replaced by the modern 
notation (r) in order to highlight that this is a Fukui 
function of second order. Its physical meaning is to reveal 
nucleophilic and electrophilic sites on a molecular system 
at the same time. Mathematically it is defined in terms of 
the derivative of the Fukui function, /(r) [14], with respect 
to the number of electrons, N. Through a Maxwell rela- 
tion, this LRD may be interpreted as the variation of rj 
(the molecular hardness which measures the resistance 
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Figure 1 A sketch of the molecular structure of the naringenin 
flavonoid. This figure shows a sketch of the molecular structure of 
the naringenin flavonoid. 



to charge transfer [21]) with respect to u(r), the external 
potential. The definition of/^(r) is shown as indicated 
by Morell et al. [5,9]: 

\dN J V(I) [Sv(t)\ n 

As mentioned above, DD allows one to obtain simul- 
taneously the preferably sites for nucleophilic attacks 
(f( 2 \r) > 0) and the preferably sites for electrophilic 
attacks (f^ 2 \r) < 0) into the system at point r. DD has 
demonstrated to be a robust tool to predict specific sites 
of nucleophilic and electrophilic attacks in a much more 
efficient way than the Fukui function by itself because dual 
descriptor is able to distinguish those sites of true nucle- 
ophilic and electrophilic behavior, in consequence some 
works have been published with the aim of remarking the 
powerfulness oif^ (r) and all those LRDs depending on 
DD [5,9,11,16-19]. 

The general working equation to obtain DD is given 
by the difference between nucleophilic and electrophilic 
Fukui function [11]. A well-known first level of approxi- 
mation implies the use of finite difference method where 
to the sum of electronic densities of the system with 
one more electron and one less electron is subtracted by 
the double of the total electronic density of the original 
system. Since this level of approximation implies a time- 
demanding computing, a second level of approximation 
has been used for some years where the densities of FMOs 
provide an easier-to-compute working equation: 

/ (2) (r) = /+ (r) -/" (r) ~ p L (r) - p H (r), (2) 

where densities of LUMO and HOMO are represented by 
p L (r) and p H (r), respectively. 

Hence, when an interaction between two species is well 
described through the use of this LRD, it is said the 
reaction is controlled by frontier molecular orbitals (or 
frontier-controlled) under the assumption that remaining 
molecular orbitals do not participate during the reaction. 

The dual descriptor can also be condensed through an 
appropriate integration within the k th - atomic domain Q^: 

f f< 2 Hr)dr=ti 2) . (3) 

jQ k 

When/^ > 0 the process is driven by a nucleophilic 
attack on atom k and then that atom acts an electrophilic 
species; conversely, when^; ; < 0 the process is driven by 
an electrophilic attack over atom k and therefore atom k 
acts as a nucleophilic species. 

Settings and computational methods 

All computational studies were performed with the 
Gaussian 09 [22] series of programs with density func- 



Glossman-Mitnik Chemistry Central Journal 201 3, 7:1 55 
http://journal.chemistrycentral.eom/content/7/1/155 



Page 3 of 7 




tional methods as implemented in the computational 
package. The equilibrium geometries of the molecules 
were determined by means of the gradient technique. 
The force constants and vibrational frequencies were 
determined by computing analytical frequencies on the 
stationary points obtained after the optimization to 
check if there were true minima. The basis set used 
in this work was MIDIY, which is the same basis set 
as MIDI! with a polarization function added to the 
hydrogen atoms. The MIDI! basis is a small double-zeta 
basis with polarization functions on N-F, Si-Cl, Br, and 
I [23-28]. 

For the calculation of the molecular structure and 
properties of the studied system, we have chosen the 
hybrid meta-GGA density functionals M06, M06L, M06- 
2X and M06HF [29], which consistently provide satis- 
factory results for several structural and thermodynamic 
properties [29-31]. All the calculations were performed 
in the presence of water as a solvent, by doing IEFPCM 
computations according to the SMD solvation model 
[32]. 

Within the conceptual framework of DFT [14,21], 
the chemical potential /x, which measures the escap- 
ing tendency of electron from equilibrium, is defined 



as: 



dE 

dN 



-x 



(4) 



v(r) 



where x is the electronegativity. 

The global hardness rj can be seen as the resistance to 
charge transfer: 



_ 1 /8 2 £\ 
71 ~ 2 \dNz) v Cr) 



(5) 



Using a finite difference approximation and Koopmans' 
theorem [25-28], the above expressions can be written 

as: 



'-(I + A) * -(€ L + e H ) 



^«i(/-A)»i(€£-€ if ) 



(6) 



(7) 



where €h and €i are the energies of the highest occupied 
and the lowest unoccupied molecular orbitals, HOMO 
and LUMO, respectively. However, within the context 
of density functional theory, the above inequalities are 
justified in light of the work of Perdew and Levy [33], 
where they commented on the significance of the high- 
est occupied Kohn-Sham eigenvalue, and proved the 
ionization potential theorems for the exact Kohn-Sham 
density functional theory of a many-electron system. 
In addition the use of the energies of frontier molec- 
ular orbitals as an approximation to obtain / and A is 
supported by the Janaks Theorem [34]. In particular, 
The negative of Hartree-Fock and Kohn-Sham HOMO 
orbital has been found to define upper and lower lim- 
its, respectively, for the experimental values of the first 
ionization potential [35] thus validating the use of ener- 
gies of Kohn-Sham frontier molecular orbital to cal- 
culate reactivity descriptors coming from Conceptual 
DFT. 

The electrophilicity index co represents the stabilization 
energy of the systems when it gets saturated by electrons 
coming from the surrounding: 



(/ + A) 2 _ (e L + e H f 



2T] 4(1 -A) 4(e L -e H ) 



(8) 
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The electron donating (co ) and electron accepting (&> + ) 
powers have been defined as [36]: 



co = 



(37 + A) 2 
16(1 -A) 



and 



(/ + 3A) 2 
16(1 -A) 



(9) 



(10) 



It follows that a larger o; + value corresponds to a better 
capability of accepting charge, whereas a smaller value of 
co~ value of a system makes it a better electron donor. In 
order to compare co + with -co~ , the following definition of 
net electrophilicity has been proposed [37]: 



A&r^ = co + — (-co ) = co + + co 



(id 



that is, the electron accepting power relative to the elec- 
tron donating power. 

Results and discussion 

The molecular structure of Naringenin was pre-optimized 
by starting with the readily available PDB structure, 
and finding the most stable conformer by means of the 
Conformers module of Materials Studio through a ran- 
dom sampling with molecular mechanics techniques and 
a consideration of all the torsional angles. The structure 
of the resulting conformer was then optimized with the 
M06, M06L, M06-2X and M06-HF density functionals in 
conjunction with the MIDIY basis set. 

The validity of the Koopmans' theorem within the 
DFT approximation is controversial. However, it has been 
shown [35] that although the KS orbitals may differ in 
shape and energy from the HF orbitals, the combination 
of them produces Conceptual DFT reactivity descriptors 
that correlate quite well with the reactivity descriptors 
obtained through Hartree-Fock calculations. Thus, it is 
worth to calculate the electronegativity, global hardness 
and global electrophilicity for the studied systems using 
both approximations in order to verify the quality of the 
procedures. 

The HOMO and LUMO orbital energies (in eV), ion- 
ization potentials I and electron affinities A (in eV), and 
global electronegativity x> total hardness rj, and global 
electrophilicity co of the Naringenin molecule calculated 
with the M06, M06L, M06-2X and M06-HF density func- 
tionals and the MIDIY basis set are presented in Table 1. 
The upper part of the table shows the results derived 
assuming the validity of Koopmans' theorem and the 
lower part shows the results derived from the calcu- 
lated vertical I and A. As can be seen from Table 1, 
the Koopmans theorem holds approximately for the den- 



Table 1 HOMO and LUMO orbital energies (in eV), 
ionization potentials I and electron affinities A (in eV), and 
global electronegativity x, total hardness rj, and global 
electrophilicity co of Naringenin calculated with the M06, 
M06L, M06-2X and M06-HF density functionals and the 
MIDIY basis set 



Property 


M06 


M06L 


M06-2X 


M06-HF 


HOMO 


-5.9155 


-4.7242 


-7.3281 


-9.2653 


LUMO 


-0.6052 


-1.4060 


0.1298 


1.4321 


X 


3.2604 


3.0651 


3.5992 


3.9166 


rj 


2.6552 


1.6591 


3.7290 


5.3487 


00 


2.0018 


2.8313 


1.7370 


1.4339 


1 


7.4000 


6.9619 


8.1214 


8.8799 


A 


0.9393 


0.7328 


0.9528 


0.8997 


X 


4.1697 


3.8474 


4.5371 


4.8898 


r\ 


3.2304 


3.1146 


3.5843 


3.9901 


00 


2.6911 


2.3763 


2.8716 


2.9962 



The upper part of the table shows the results derived assuming the validity of 
Koopmans' theorem and the lower part shows the results derived from the 
calculated vertical I and A. 



sity functionals which include some percentage of HF 
exchange, but it fails in part for the M06L density func- 
tional (without inclusion of HF exchange) (for the elec- 
tronegativiy and total hardness but not for the global 
electrophilicity). 

The condensed Fukui functions can also be employed 
to determine the reactivity of each atom in the molecule. 
The corresponding condensed functions are given by 
f+ = q^(N + 1) — qk(N) (for nucleophilic attack),^, - = 
qk(N) — qk(N — 1) (for electrophilic attack), and f® = 
[ qk(N + 1) — qk(N — 1)] /2 (for radical attack), where q% 
is the gross charge of atom k in the molecule. 

It is possible to evaluate condensed Fukui functions 
from single-points calculations directly, without resorting 
to additional calculations involving the systems with N-l 
and N+l electrons: 



aek 



c ai H~ c ai ^ ] c bi&ab 



and 

/*- = £ 

aek 



c ai / , c bi^ab 
b^a 



(where i = LUMO) 
(12) 

(where i = HOMO) 
(13) 



with c a i being the LCAO coefficients and S a b the overlap 
matrix. The condensed Fukui functions are normalized, 
thus Ekfk = 1 and/? =[/+ +/,"] 12. 
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The condensed Fukui functions have been calculated 
using the AOMix molecular analysis program [38,39] 
starting from single-point energy calculations. We have 
presented, discussed and successfully applied the des- 
cribed procedure in our previous studies on different 
molecular systems [40-43] . 

The condensed dual descriptor has been defined as 
f( 2 \r)k = fk ~fk [5>9]. From the interpretation given to 
the Fukui function, one can note that the sign of the dual 
descriptor is very important to characterize the reactiv- 
ity of a site within a molecule toward a nucleophilic or an 
electrophilic attack. That is, if/^(r)^ > 0, then the site is 
favored for a nucleophilic attack, whereas if/^(r)^ < 0, 
then the site may be favored for an electrophilic attack 
[5,9,44]. 

The electrophilic f~ and nucleophilic f + condensed 
Fukui functions and/^ 2 ^(r) over the atoms of the Narin- 
genin molecule calculated with the M06, M06L, M06-2X 
and M06-HF density functionals and the MIDIY basis 
set are shown in Table 2. The actual values have been 
multiplied by 100 for an easier comparison. 

It can be concluded from the analysis of the results 
on Table 2 that the M06, M06L, M06-2X and M06-HF 



density functionals predict that CIO will be the preferred 
site for nucleophilic attack. The four density functionals 
considered in this study display a large negative value of 
the condensed dual descriptor f^ 2 \v) over 018, imply- 
ing that this will be the preferred site for the electrophilic 
attack. 

The electrodonating (co~) and electroaccepting (&> + ) 
powers and net electrophilicity Ao; ± of the Narin- 
genin molecule calculated with the M06, M06L, M06- 
2X and M06-HF density functionals and the MIDIY 
basis set are presented in Table 3. The upper part 
of the table shows the results derived assuming the 
validity of Koopmans' theorem and the lower part 
shows the results derived from the calculated vertical I 
and A. 

The results from Table 3 clearly indicate that Narin- 
genin is an electrodonating molecule, with the same result 
predicted by all the four density functionals considered 
in this study. However, although the tendency is the 
same, the results for these descriptors are in poor agree- 
ment between those calculated assuming the validity of 
the Koopmans' theorem, and those coming from energy 
differences. 



Table 2 Electrophilic f and nucleophilic f + condensed Fukui functions and f (2) (r) over the atoms of the Naringenin 
molecule calculated with the M06, M06L, M06-2X and M06-HF density functionals and the MIDIY basis set 







M06 






M06L 






M06-2X 






M06-HF 




Atom 


f" 


f" 




f+ 


f" 




f+ 


f" 




f+ 


f" 




1 0 


2.13 


0.99 


1.14 


2.32 


0.35 


1.97 


1.58 


1.23 


0.35 


1.05 


0.50 


0.55 


2C 


11.39 


0.30 


11.09 


10.86 


0.29 


10.57 


11.68 


0.28 


11.48 


11.58 


0.06 


11.52 


3C 


1.36 


1.26 


0.10 


1.27 


1.21 


0.06 


1.45 


0.99 


0.46 


1.62 


0.88 


0.74 


4C 


0.61 


0.73 


-0.12 


0.51 


0.32 


0.19 


0.82 


1.14 


-0.32 


1.38 


0.23 


1.15 


5C 


6.30 


6.30 


0.00 


5.23 


8.06 


-2.83 


8.31 


0.91 


7.40 


11.67 


0.27 


11.40 


6C 


0.36 


9.08 


-8.72 


0.30 


10.14 


-9.84 


0.34 


3.45 


-3.11 


0.33 


2.46 


-2.13 


7C 


0.49 


8.61 


-8.12 


0.41 


2.57 


-2.16 


0.51 


2.67 


-2.16 


0.63 


2.75 


-2.12 


9C 


17.38 


-0.01 


17.39 


16.22 


-0.01 


16.23 


18.71 


0.03 


18.68 


16.60 


0.03 


16.57 


10C 


21.29 


1.27 


20.02 


22.36 


1.07 


21.29 


19.76 


0.87 


18.89 


20.73 


0.78 


19.95 


11 C 


9.32 


0.21 


9.11 


9.57 


0.15 


9.42 


9.10 


0.31 


8.79 


9.04 


0.09 


8.95 


14C 


0.05 


2.43 


-2.38 


0.07 


1.05 


-1.08 


0.06 


5.41 


-5.35 


0.07 


6.81 


-6.74 


15C 


0.59 


1.70 


-1.11 


0.42 


0.81 


-0.39 


0.55 


3.29 


-2.74 


0.64 


4.18 


-3.54 


160 


4.11 


0.02 


4.09 


4.44 


0.03 


4.41 


3.45 


0.02 


3.43 


2.68 


0.01 


2.67 


17C 


1.02 


1.03 


-0.01 


0.69 


0.54 


0.15 


1.53 


0.84 


0.69 


2.24 


0.11 


2.13 


180 


18.58 


44.09 


-25.51 


20.18 


65.26 


-45.08 


17.95 


45.56 


-27.61 


16.24 


41.45 


-25.21 


190 


2.31 


0.98 


1.33 


2.65 


1.10 


1.55 


1.81 


0.40 


1.41 


1.28 


0.06 


1.22 


20 C 


0.20 


3.63 


-3.43 


0.10 


0.90 


-0.80 


0.19 


9.54 


-9.35 


0.26 


10.09 


-9.83 


23 C 


0.05 


3.72 


-3.67 


0.02 


0.91 


-0.89 


0.03 


10.56 


-10.53 


0.03 


10.99 


-10.96 


28 C 


0.32 


5.50 


-5.18 


0.26 


1.73 


-1.47 


0.34 


15.69 


-15.35 


0.42 


20.75 


-20.33 


30 0 


0.08 


7.15 


-7.07 


0.08 


2.19 


-2.11 


0.06 


15.62 


-15.56 


0.05 


12.24 


-12.19 



The actual values have been multiplied by 1 00 for an easier comparison. H atoms are not shown. 
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Table 3 Electrodonating (co ) and electroaccepting (co + ) 
powers and net electrophilicity Aw ± of Naringenin 
calculated with the M06, M06L, M06-2X and M06-HF 
density functionals and the MIDIY basis set 



Property 


M06 


M06L 


M06-2X 


M06-HF 


0J~ 


3.9637 


4.5712 


4.0026 


4.0609 


co + 


0.7035 


1.5061 


0.4035 


0.1443 




4.6672 


6.0773 


4.4061 


4.2052 


0J~ 


6.2176 


4.6892 


5.5882 


5.9398 


C0+ 


1.2124 


0.8419 


1.0511 


1.0500 




7.4300 


5.5311 


6.6393 


6.9898 



The upper part of the table shows the results derived assuming the validity of 
Koopmans' theorem and the lower part shows the results derived from the 
calculated vertical I and A. 



Conclusions 

From the whole of the results presented in this contri- 
bution it has been clearly demonstrated that the sites of 
interaction of the Naringenin molecule can be predicted 
by using DFT-based reactivity descriptors such as the 
hardness, softness, and electrophilicity, as well as Fukui 
function calculations. These descriptors were used in the 
characterization and successfully description of the pre- 
ferred reactive sites and provide a firm explanation for the 
reactivity of the Naringenin molecule. 

The M06 family of density functionals (M06, M06L, 
M06-2X and M06-HF) used in the present work leads to 
the same qualitatively and quantitatively similar descrip- 
tion of the chemistry and reactivity of the Naringenin 
molecule, yielding reasonable results. However, for the 
case of the M06-2X and M06-HF density functionals, 
which include a large portion of HF exchange, the calcu- 
lations considering the validity of the Koopmans' theorem 
lead to negative electron affinities. 

The calculated descriptors are in agreement with the 
known experimental facts about the chemical reactivity 
of the Naringenin molecule presented in the literature 
(with the exceptions mentioned on the paragraph above). 
Thus, this make us confidents that similar studies can 
be pursued with the same degree of accuracy on another 
flavonoids with analogue structures. 
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